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G\ ■ Abstract 

OV 

i—* \ We consider deuteron formation in heavy ion collisions at intermediate 

^O ■ energies. The elementary reaction rates (Nd — ► NNN etc.) in this 

r^r . context are calculated using rigorous Faddeev methods. To this end an 

in-medium Faddeev equation that consistently includes the energy shift 
and Pauli blocking effects has been derived and solved numerically. As a 

t— i ■ first application we have calculated the life-time of deuteron fluctuations 

for nuclear densities and temperatures typical for the final stage of heavy 
ion collisions. We find substantial differences between using the isolated 

(^ ■ and the in- medium rates. 

(N' 

o: 

9£ ' 1 Introduction 

On ; 

,-G ■ Nuclear matter is an example of a strongly correlated many particle system. 

One prominent consequence is the formation of bound states (clusters, frag- 
ments) observed in heavy ion reactions. Here we address the formation of 
deuterons at intermediate energies, i.e. for E/A < 200 MeV/u. Within the 
quantum statistical approach to describe the complicated dynamics we employ 
the Green function method [lj. The cluster mean field approximation de- 
couples the hierarchy and leads to rigourous few-body equations for the two-, 
three-, four-particle correlations. This method is tiedly connected to the self 
consistent RPA approach extended to finite temperatures ||. 

Deuteron formation is directly related to the Nd — > NNN break-up cross 
section. Photoinduced reactions have also been considered [[§]. Because of the 
energies considered, pion induced reactions can be neglected. 

Treating deuteron formation within the cluster Hartree-Fock approxima- 
tion allows us to consistently include all medium modifications as they appear 
in the respective two- and three-body equations. These are the self energy and 
the Pauli blocking effects. 
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2 Reactions 

The quantity of interest in the quantum statistical approach is the general- 
ized quantum Boltzmann equation for the nucleon f^, deuteron fa, etc. dis- 
tributions. Here, we consider the collision integral to show the relevance of 
three-body reaction rates. The feeding of the nucleon density is driven by the 
collision integral (see e.g. Ref. |5] for an application to heavy ion collisions) 

JN(p,t) = JZ Ut (p,t)f N (p,t)-J™(p,t)f N (p,t), (1) 

where we have used /at = 1 — /jv- To be more explicit we give J7^ ut (p, t), 

J^(p,t) = J d 3 k J d'hd'hKkplTlhh)] 2 f N (h,t)f N (k 2 ,t)f N (k,t) 

+ f d 3 k I d 3 k t d 3 k 2 d 3 k 3 {(kplUolhhk^l 2 

x/jvOi, t)f N (k 2 , t)f N (k 3 , t)f d (k, t) 
+ ... (2) 

where dots stand for other possible contributions mentioned before. The quan- 
tity Uq appearing in (|2]) is the break-up transition operator for Nd — > NNN. 
For the isolated three-body problem Uq determines the break-up cross section 
cr via 

a = , ~ fd 3 k 1 d 3 k 2 d 3 k 3 \(kp\U \k 1 k 2 k 3 )\ 2 

x2ir5(E f - Ei) (2n) 3 5 {3) (k l + k 2 + k 3 ). (3) 

So far the strategy has been to implement the experimental cross section into 
the above equation. This has then been solved, for a specific heavy ion colli- 
sion 0. Using experimental cross sections respectively isolated cross sections 
may not be sufficient in particular in the lower energy regime. The cross 
section itself depends on the medium, e.g. blocking of internal lines or self 
energy corrections of the respective three-body Green functions. To this end 
we have derived a three-body Faddeev type equation || 0]. We use the AGS 
formalism || for the three-body algebra and solve the respective equations nu- 
merically. The equation for the three-particle Green function derived within 
the cluster Hartree-Fock approximation reads 

r f v AM + fihh , (/i/2-/i/2)V(12) + cycl.perm. 

Z — E\ — E 2 — £ 3 Z — E\ — £ 2 — E 3 



where e denotes the quasi particle energies evaluated in Hartree-Fock approx- 
imation and /i = f(ei) the Fermi functions f(e) = (exp[/3(e — jj) + l]~ x with 
the inverse temperature (3 and the chemical potential fi (for the time being 
we assume symmetric nuclear matter). Here, we use equilibrium distributions 
to solve Faddeev equations. This is justified within the linear response the- 
ory, where nonequilibrium quantities are expressed through equilibrium ones, 
because of small fluctuations only. The question of self consistency has been 
addressed for the much simpler case of two-particle correlations, e.g. in Ref. ||. 
Within the AGS formalism (extended here to finite temperatures and den- 
sities) the break-up operator Uq is simply related to the elastic/rearrangement 
scattering amplitude U a p connecting the channels^ j3 — > a. It is therefore 
sufficient to present the AGS equation for the transition operator U a p only, 
viz. 



U a p(z) = 5 a/3 



/1/2/3 + /1/2/3 

Z — E\ — 62 — £3 



+ J2 T ^ hhh + hhh U lP {z) (5) 



-,{j) J1J2J3 -r J1J2J3 



with 5 = 1 — 5 and / = 1 — /. The two-body t matrix T 3 has been solved 
on the same footing consistently including all medium effects, i.e. T 3 is the 
solution of the in-medium two-body problem, e.g. for 7 = 3 

T 3 (3) = (1 - f 3 + 9(6, + e,))- 1 V 2 + V 2 hh ~ hh T 3 (3) , (6) 

z — S\ — 82 — £3 

where g(uf) = (exp[(3(uj — 2/x) — 1] _1 is the Bose function for two nucleons. 

3 Results 

The AGS equations have been solved for a separable Yamaguchi potential. 
To get an impression of the quality of the calculation the isolated cross sec- 
tion is given in Fig. [3] along with the experimental data on neutron deuteron 
scattering JRJ. 

From inspection of Fig. |2] we see that the in-medium cross section is sig- 
nificantly enhanced compared to the isolated on. The threshold is shifted 
to smaller energies, which is because the binding energy of the deuteron be- 
comes smaller. We observe that for higher energies the medium dependence 
of the cross section becomes much weaker, which a posteriori justifies the use 



■l- The channel notation a, (3 = 1,2,3 labels the respective spectator in the three-body 
system. 
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Figure 1: A comparison of the to- 
tal, elastic, and break-up cross sec- 
tions nd — > nd, nd —>■ nnp with the 
experimental data of Ref. |TD| . 



Figure 2: In-medium break-up cross 
section at T = 10 MeV. Isolated cross 
section is shown as solid line, other 
lines show different nuclear densities. 



of isolated cross sections (along with the impulse approximation) when higher 
energies are considered ||. 

From linearizing the Boltzmann equation it is possible to define a break-up 
time for small fluctuations of the deuteron distributions. For small fluctuations 
Sf(t) = fd(t) — f% from the equilibrium distribution f% linear response leads 
to 

d t 6f d (P,t) = —6f d (P,t) (7) 

Tbu 

where the "life time" of deuteron fluctuations has been introduced , 

4 



r bu 



-J dk N d 3 hd 3 k 2 d 3 k 3 \(kp\U \kik 2 k 3 )\ 2 fJ 2 f 3 f(k N )27r8(E-E ). (8) 

which can be related to the break-up cross section given in Fig. || For low 
densities the life time (as a function of the deuteron momentum P) and the 
inverse life time, i.e. the width, at P = along with the deuteron binding 
energy for comparison is shown in Fig. |3|. These times have to be compared to 
the approximate duration of the heavy ion collision of about 200 fm. 

Another important time scale is the chemical relaxation time for small 
fluctuations of the deuteron density 8nd(f) = nd{t) — n d from the equilibrium 
distribution n d . Using detailed balance and linearized rate equations the re- 
laxation time is given through 

d 



j.Sn d (t) = 5n d (t). 

dt T re l 



(9) 
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Figure 3: Deuteron break-up time at 
T = 10 MeV and nuclear density n = 
0.007 fm~ 3 . Solid line with medium 
dependent cross section a given in 
Fig. H Short dashed with P cm = 
and dashed line isolated cross section. 



Figure 4: Relaxation time for small 
fluctuations of the deuteron density 
from chemical equilibrium at a tem- 
perature of T = 10 MeV. Line coding 
as in Fig. |] 



The basic quantity driving the time scale is again the break-up cross section 



Trcl 



— = / d 3 p N d 3 p d f(k N )g(p d )\v d - v N \ a (E) 



n% + An , 



o ^o 



n N n d 



(10) 



The resulting relaxation time as a function of the uncorrelated nuclear density 
is given in Fig. |j. 



4 Conclusion and Outlook 

Our results show that medium dependent cross sections in the respective colli- 
sion integrals lead to shorter reaction time scales. Chemical processes become 
faster. This also effects the elastic rates that are related to thermal equilibra- 
tion. 

The basis of this result is the cluster Hartree-Fock approach that in our 
approximation includes correlations up to three particles in a consistent way. 
The equations driving the correlations are rigorous. The respective one-, two- 
and three-body equations are solved, in particular for the three-particle case 
Faddeev/AGS type equations have been derived in Ref. || |7|]. The AGS ap- 
proach is particularly appealing since it allows generalizations to n-particle 
equations in a straight forward way. Results for the three-body bound states 



in medium will be published elsewhere [fTI]. As expected form the deuteron 
case, the triton binding energy changes with increasing density up to the Mott 
density, where E t — 0. 

The production rates, spectra etc. of light charged particles in heavy ion 
collisions at intermediate energies may change because of the much smaller 
time scales induced through the medium dependence compared to the use of 
free cross sections (respectively experimental cross sections). To this end some 
notion of the relevant densities and temperatures during the heavy ion collision 
(during the final stage) should be achieved. 
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